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ABSTRACT 

The 7r°77 vertex for virtual photons of squared masses q\ and q\ plays a vital role 
in several physical processes; for example for qf < 0, q\ < 0, in the two-photon physics 
reaction e + e~ —>■ e + e~7c°, and for q\ > 0, q\ > 0, in the annihilation process e + e~ — > 
7r°l + l~ . It is also of interest because of its link to the axial anomaly. We suggest a new 
approach to this problem. We have obtained a closed analytic expression for the vertex in 
the limit in which at least one of | | and \q%\ is large for arbitrary fixed values of the ratio 
q\jq\- We compare our results with those obtained previously by Brodsky and Lepage. It 
should be straightforward to test our predictions experimentally. 
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1. Introduction 



As is well known the famous Adler-Bell-Jackiw anomaly leads to a non- vanishing am- 
plitude for neutral pion decay, n° — > 27, in the chiral limit when the pion mass and the 
current masses of light quarks go to zero: m n — > 0, m M ,m,j — > 0. The result is a conse- 
quence of the linear divergence of the triangle diagram and the whole contribution actually 
comes from short distances when a gauge invariant regularization is utilized. A natural 
question arises: does anything change when the photons are off-shell? At first glance, since 
the anomalous contribution has a point-like character, nothing changes seriously provided 
the virtualities (masses) of photons remain smaller than the scale of regularization, which, 
however, can be chosen arbitrarily large. This simple observation led Jacob and Wu some 
years ago to very strong statements about certain physical processes. For instance, if 
the decay amplitude Z — > n° + 7, which is intimately related to n° —>■ 7* + 7, is controlled 
by a point-like amplitude and does not decrease as 771(7*) — > m z, one has a 7r°7 production 
rate wildly incompatible with LEP data. The above expectation seems, of course, quite 
unintuitive per se; still stranger that at the time of publication of ref. the correct an- 
swer had been known for some years. In fact an extremely small value of the amplitude 
of Z° — > 7r° + 7 decay obtains because of the decrease of the amplitude tt° — ■> 7 + 7* as 
the virtuality of the photon increases |§ @] H- Another similar amplitude, which had 
also been considered earlier, is n° — > 7* + 7* with equal, large negative virtualities for the 
photons, q\ = q\ = —Q 2 @. For this case a simple operator product expansion (OPE) 
leads to the natural conclusion that the amplitude falls as 1/Q 2 in the asymptotic region 
Q 2 -> 00. 

The paper |J has already been criticized 0] but since to some extent it triggered 
this investigation we shall give in the next section a simple explanation of the "paradox" 
described above. 

There are interesting experimental motivations to study the amplitude n° — > 7* + 7* 
for non-vanishing masses q 2 7^ 0, q\ 7^ of the virtual photons. This amplitude can 
be measured in principle as a function of negative q\ < and q 2 < in "two-photon 
physics" (see Fig.l). For positive values qf, q\ > it enters the amplitude of the process 
e+ + e~ -> 7t° + /+ + I- (Fig.2). Though this process is too small to be measured when it 
goes through the Z boson it can be observed in e + e~ collisions at lower energies. 

Consequently the interest in the amplitude 7r° — » 7* + 7* is far from exhausted and 
there is a growing literature on the subject || , || , [|10| . Indeed experiments in "two- photon 



physics" have now been started by CLEO at Cornell. For a summary of data up to the 
end of 1994 see ref. 0. 

The Feynman amplitude for % (k) — > j*(qi) + 7*(— 92) can be written in the form 

M = M(ql q\) ■ e^ XtT E»{ qi y{-q 2 )q X k° (1.1) 

where 

q={qi + q 2 )/2, k = qi -q 2 , and q 2 = (q 2 + q 2 )/2 - m 2 jA. (1.2) 
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Two results, of great importance, have been known for many years. In the case of equal 
space-like momenta, q\ = q\ < the operator product expansion (OPE) yields the follow- 
ing behaviour as q\ = q\ — > — oo 

M(gV) = M(0,0)(--^J for g 2 ^-oc (1.3) 



where we have used q 2 ~ qf — q\ from eqs.( |1.2|) and where f n ~ 93MeV. The amplitude 
M(0, 0) is the usual amplitude for it — > 7 + 7 into real photons: 

M(0,0)=^.i. (1.4) 

For the case of one photon on mass-shell (q 2 = 0), and the other with large positive 
virtuality (q 2 > 0), Brodsky and Lepage have used the QCD infinite momentum frame 
approach for the description of the pion wave function and have obtained the asymptotic 
behaviour 

M(2g 2 ,0) = M(0,0) ( — for q 2 -> 00 (1.5) 

where we have used g 2 ~ ?i/2 from eqs.(|L~2|). 

Generally the amplitude M is a function of the two variables q\ and q\. It is interesting 
to consider the case when q\ and q\ both go to infinity, but their ratio remains finite. As 
far as we know for this case no reliable results exist. 

In this paper we try to derive some conclusions about the asymptotics in the latter 
region, attacking the problem by means of a triangle diagram in which the pion-quark- 
antiquark vertex is given in its most general (non-pointlike) form. In the next section 
we shall argue that this diagram, though mixing hadron and quark degrees of freedom, 
gives reliable results for our purposes. This is not a trivial point and the straightforward 
analogy with the triangle diagram always used for 7r° decay into real photons is somewhat 
misleading. In the latter case one actually does not use a diagram containing the pion and 
the quarks, but rather a diagram containing the axial current. It is only through PCAC 
that we can connect it to n° — > 27 decay. The axial current vertex is pointlike by definition, 
quite different from the case of the pion vertex which is considered below. 

Our results are the following. First we have tested our approach by reproducing the 
OPE asymptotics for the case of large equal and negative virtualities, q\ = q 2 < 0, eqs. 

Second, for the more general case of a fixed ratio q\jq\, but q\, q\ going to infinity we 
find a "scaling law" : 

8^' 



M(q(,qi) = M(0,0) ( -- ^ ) </>( w ), (1. 6 ) 



2 2 
_ Ti ~ Q2 

Ul —9 7y 

QI + Q2 
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where 0(w) is a certain dimensionless function. 

Note that the factor f 2 /q 2 appears naturally in the particular cases described by 
eqs. ( |1.3|) and ( |1.5| ). The variable ui lies in the interval — 1 < ui < 1 since we consider the 
physically interesting cases of either both q 2 's negative, q\,q\ < 0, (Fig.l) or both positive, 
q\,ql > (Fig.2). Eqs. flL3|) and (|L5|) imply that 

0(0) = 1, 0(1) =3/2. (1.7) 

The principle aim of this paper will be to calculate 0(w) for all — 1 < u < 1. We shall be 
able to find 0(w) in the chiral limit, — > 0, and under some natural physical assumptions. 
Our result for 0(w) reads 

= — [(1 + £j)ln(l +u) + {l- u)ln(l - u)]. (1.8) 

At u = 0(0) = 1, as required by eq.(|L7|); at u = 1 0(1) = 2 In 2 « 1.4. The latter value 
is slightly different from that given by eq. ( |1.7|) which is not surprising since our approach 
differs from that of ref. ||. 

In fact the form ( |1.6| ) is not new (see, for example, 0). Generally speaking the 
exact form of the function 0(w) is determined by the detailed form of the pion wave 
function. Our calculation shows, however, that one can obtain an approximate expression 
for 0(w) without an explicit knowledge of the pion wave function, provided one assumes 
reasonable asymptotic properties of the wave function or the related form factors. Actually 
in our computation we use an expansion with respect to a formal parameter related to the 
asymptotic behaviour of pion wave function. The leading behaviour of the form factor 
with respect to this formal parameter turns out to be independent on details of the pion 
wave function. In ref. [[| the validity of the treatment is restricted to the region of small 
< 0.5 (for a further analysis in the sum rule approach see [1(|). I n re f- flUt a series for 



4>{oj) is obtained in which succeeding terms are suppressed by powers of log Q 2 but it is 
claimed that the result is only reliable for small omega. Contrary to the approaches of refs. 



|J and [[12] we use a numerical approximation, which, we hope, may be valid for |u;| < 1. 
Note also that our treatment ignores corrections of order l/logQ 2 - Therefore our result 
may be formally interpreted as an expression for the amplitude in the limit of very large 
Q when the corrections ~ 1/logQ 2 can be neglected. However since we exploit a quite 
different physics as compared to ref. 0] we think that our result is a good approximation 
also for 'intermediate' Q which are larger than typical hadronic scale but logQ 2 is not 
negligible. 

In the last part of the paper we apply our formalism to the it decay into real photons. 
We calculate the decay amplitude and compare it with the PCAC result. This leads to 
a relation between one of the form factors describing the pion-quark-antiquark vertex at 
small momenta and /„-. The implication of this relation is that at small momenta the 
effective pion-quark-antiquark Lagrangian has the form 

L eff = ^(d^)(qr 3 ^ 5 q). (1.9) 
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The paper is organized as follows. In section 2 we discuss the Jacob- Wu "paradox" and 
introduce the main formalism which we use. In section 3 we calculate the constant f n in 
terms of integrals over the introduced form factors and study the asymptotics of n° — > 
7* + 7* at equal large photon masses qf = q\ = q 2 — > ±oo. The OPE result for this 
asymptotics is reproduced. In section 4 we calculate the function (f)(u>), eq. (|1.8|) . In 
section 5 we deal with 7r° decay into real photons and derive the relation between the zero 
momentum value of one of our form factors and f n , corresponding to eq. (|1.9|) . In the 
conclusion we summarize and discuss briefly our results. 



2. General formalism 

To set the stage we begin by discussing the Jacob- Wu "paradox" in which the tv° — > 
7* + 7* amplitude does not decrease with the large photon virtualities. 

In considering it decay into two real photons one uses the PCAC relation and the 
divergence of the axial- vector current A M , corrected by the Adler-Bell-Jackiw anomalous 
term 

2 _ 

d^A^Uml^ + ^F^F^, A^ = q^ W q (2.1) 

where the notation is standard. 

From eq.(|2.1|) one gets, by well-known arguments, for the n — > 27 amplitude 



M = < 2 7 |^|0 > +±- • -f^ ■ e^ x ^{ qi y{-q 2 )q X k°. (2.2) 

The matrix element of d^A^ between the vacuum and 27 state is proportional to k 2 where 
k is the four momentum of the pion and of the 27 state. Since there are no relevant 
physical massless particles in the limit k^ — > the first term in eq.( |2.2j ) vanishes and we 
remain with the second term which is the famous result for the tt — > 27 amplitude. Why 
is this not correct when the 7s are off shell? After all, as has already been mentioned, 
the contribution proportional to F^F^ comes from very short distances and therefore 
can not change unless the masses of the photons reach the regularization scale, which is 
arbitrarily large. 

To understand the changes appearing when the virtualities of the photons are not zero, 
q\ ^ 0, q\ 7^ 0, let us go back one step and write down the full expression for < 27|<9 M A^|0 > 
in eq. (J2.2[), for light but not completely massless quarks, m q = m u = rrid 7^ 0. One has W3 



< 2 7 |<V^|0 >= -— , ■ e^ x ^{ qi y{-q 2 )q X k°* 



(2.3) 



1 + 2m g dx dy 
'0 Jo 



l-x 



(qfx + q 2 y)(l — x — y) + xym\ — m 2 



The first term in this equation (unity in the parenthesis) corresponds to the short distance 
contribution. It can be derived, say, by introduction of a Pauli-Villars regulator fermion. 
After regularization it is permissible to use the equations of motion for the operators in 
calculating the divergence of A^. One then obtains the second term in eq.( [2.3| ) from the 
convergent triangle diagram with light quarks and with a pseudoscalar vertex 2im q r y^ l . 

Of course, one can not take expression ( |2.3j ) too seriously since the main contribution 
to the second term is determined by the small momentum domain (of order of m q ) of 
integration. This corresponds to long distances where the perturbative approach is not 
valid. However one can use eq.( |2.3| ) to resolve the paradox at least at the qualitative level. 

Suppose that due to nonperturbative long distance effects m q in eq. ( |2.3|) is a quantity 
which is closer to a constituent quark mass value rather than to the current quark mass. 
Then, in the chiral limit, we can neglect in the denominator in eq. (|2.3|) . For the case 
9i = I2 = the quark mass m q cancels out, the integral equals —1 and < 27|d M A Al |0 >= 
0. According to eq.( [2.2| ) one immediately arrives at the usual result for the tt° — > 27 
amplitude. (Note that before we simply argued that matrix element < 27|<9 M A At |0 > 
vanishes due to the absence of massless physical states.) 

For qf , q\ > m q the situation changes drastically. The second term no longer cancels 
the unit term in the brackets in eq.( |2.3|) and therefore the matrix element of d^A^ in 
eq.(|2.2|) cannot be neglected. Moreover, for qf, q\ » vn? q the integral in eq. (|2.3| ) is 
small compared to 1 and < 27*|<9 M A M |0 >^ -e 2 /4tv 2 ■ e tlu \ cr € ll e iy q x k cr . The amplitude 
A4(7T° — > 27*) goes to zero as qf, — ^ 00 as ~ 1/gf, ~ roughly speaking. 

We shall now attempt to attack the problem in a more general and rigorous fashion. 
A key object which we shall use throughout the paper is the momentum space proper 
vertex T a p(k, I) describing the transition of tt° into a light (up- or down-) quark-antiquark 
pair (Fig. 3). The formal definition of T(/c, /) is as the Fourier transform of the co-ordinate 
space vertex r(xi, x 2 ) 

i(2n) 4 5 4 (K - k) ■ r(Jfe, I) = J d 4 Xl d A x 2 e *(W2)zi-^-V2)* 2 f x 2 ) = (2.4) 

where X = (x\ + X2)/2, x = x\ — X2, and where the function r(#i, x-z) is defined as the 
matrix element of a product of Heisenberg operators of the quark fields between vacuum 
and pion states. For definiteness we use the it-quark 

f Q , )/ g(xi,x 2 ) = Waa'Wpp < K°\Tu a i {x{)up> (x 2 ) |0 > (2.5) 

where a and (3 are spinor labels; the pion state < 7r°| is taken with a definite momentum 
K. For the <i-quark the amplitude r(xi, x 2 ) is the same but of opposite sign. The equation 
( |2.5|) is written for a given colour, i.e. a summation over colour is not implied. The vertex 
T(xi,x 2 ) is related to the quark-antiquark component of the wave function of the pion, 
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but since we require a proper vertex part we have amputated the fermion legs in the 
corresponding Feynman amplitude by applying the operators and <ff 2 ^ in eq. Q2.5p . 

Using the definition (|2.5| ) and invar iance under charge conjugation one can easily show 
that for the matrix in spin space 



(2.6) 



T(x 2 ,xi)=C T (xi,x 2 )C 

where C = i'j ^ 2 is the charge conjugation matrix. 
The most general Lorentz structure of T(k, I) is 

T{k, I) = A7 5 + iS 7 A 7 5 /c A + iC(k ■ /) 7 A 7 5 / A + Da^^kxl^. (2.7) 

Here the scalar form factors A, B, C, and D depend on two scalars, I 2 and k ■ I [k 2 = m^). 

The relation ( |2.6p implies that all form factors A, B, C and D are even functions of 
/. Sometimes it is more convenient to use instead of the scalars I 2 and k ■ I the masses of 
the quark lines 



m-, 



(l + k/2) 2 , m 2 = (l-k/2) 2 . 



(2.1 



The eveness under / — > — / implies that T(k, /) is invariant under m 2 

Before going into the detailed calculations of the asymptotic behaviour of n° — ► 7* +7* 
using r(/c, /) we comment on its relationship to and on the question of validity of our 
approach. 

Using q for the isospin doublet (u, d) , f n is defined conventionally, by 



-iUk" =< 7t°(k)\q^ l5 ^-q\0> 



(2.9) 



wherein we sum over colour. This can be written in terms of our proper vertex by rein- 
stating the amputated legs. One has from eq. (|2.9|) the exact result 



if,k» = -( 7 M 75) a/3 J 



d A x 



colour ' 



1 



1 



Wi) \tK ' ( ^ l)( ^ 2) < ^M^aMlo > 



x=o 
(2.10) 

corresponding to the coordinate space representation of the diagram in Fig. 4. Thus finally 



3 / d 4 xTt 



(2.11) 



x=o 



This result does not rely on perturbation theory. The free massless quark propagators 
(l/i^)i,2 appear in equation (|2.11|) only because of our definition of T. 

For the n° — ► 7* + 7* amplitude we shall single out the diagram depicted in Fig. 5. 
Unlike the case of the calculation of f n the claim that this diagram is the most important 
one is by no means trivial. As above, the quark propagators carrying the momenta l±k/2 
appear as a result of our definition of T(k, /), so that here the use of free quark propagators 
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does not necessarily imply a perturbative approach. The free propagator corresponding 
to momentum I — q is meaningful if the quark line has a large virtuality, (I — q) 2 — > oo, 
which will be the case in our calculation. Other diagrams, for example Fig. 6 with an 
additional gluon, should be less important. The reason is that if q is a large momentum 
(or, more accurately, q 2 goes to infinity) then the biggest contribution arises when the 
flow of this momentum takes place through the minimal number of propagators. In fact 
if in Fig. 6 the upper 4-point proper Green function is a rapidly decreasing function of 
the virtualities of the external particles, the large momentum q will flow through the two 
quark propagators at the bottom of the diagram. The overall behaviour of the diagram 
will be ~ l/<? 3 or ~ l/<? 4 , instead of ~ l/<? 2 for the diagram in Fig. 5. In other words 
a soft hadronic system (the pion) will not tolerate a large momentum transfer because of 
the rapid fall of the form factors with increasing momentum transfer. 

Another class of diagrams are the radiative corrections to Fig. 5 such as shown in 
Fig. 7. In this case the flow of large momentum q does not lead to any additional power 
decrease in q 2 because of the divergence of the internal triangle diagram. The vertex 
corrections replace the QCD coupling by the running coupling evaluated at ~ q 2 . Due to 
asymptotic freedom this produces only logarithmic corrections to the diagram of Fig. 5. 

It should be stressed that the diagram of Fig. 5 can be only relied upon for large 
momentum q 2 . This is obvious when there is just one large space-time scale determined 
by q\ = q 2 = —q 2 = Q 2 . In this case the characteristic time interval between the emissions 
of the two photons is short, ~ 1/Q 2 , and no additional interactions of the quark can take 
place. This is certainly not true for small q 2 . 



3. Calculation of f n and asymptotics at large qf 
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Passing to the momentum representation in eq.( |2.11j ) and substituting the general 
structure for T(k,l) (eq. (|2.7|) ) one obtains the following expression for f n 



12 



dH 

i(2Tv) z 



B + -(k-l)C 



(I -k/2)" 
. (I - k/2) 2 



B--{k-l)C 



(l + k/2Y 
J (l + k/2)' 



• (3-1) 



The form factors A and D disappear when the trace is calculated. We change I — > — I in 
the second integral, use the fact that B and C are even functions of /, and can substitute 
Ifi — ► kfj,(k ■ l)/k 2 under the integral sign to obtain 



fir = "12 



(in 

z(2tt) z 



B+-(k-l)C 



, (I -k/2) 2 



1-2 



(fc-0 

k 2 



(3.2) 



For the amplitude 7*(?i) — > 7* (92) + ^°(k) the direct calculation of the diagram of Fig. 5 
yields the result (leaving out the polarization vectors) 
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M,t V (q, k) = 4e 2 e^xa 



dH 



i(27r) 4 (/ - q)< 



q x k a -l x k a + 2q x l 



q x k a - l x k a - 2q x l c 
(/ - k/2) 2 

B- l -{k-l)C 



B + -(k-l)C 



+ (3.3) 



(l + k/2) 2 

By changing I — > — / and g — > — g in the second integral we can represent the amplitude 
M. in the form 

M^q, k) = 4e 2 e^ X(T [A X(T (q, k) - A x °(-q, fc)], (3.4) 

where 



A ACT (g,fc) 



d 4 Z 



1 q x k a -l x k a -2q x l a 



B + \(k-l)C 



i(2n) 4 (l-q) 2 (l-k/2) 2 
The integrals of Z 17 multiplied by scalar factors must yield results proportional to q a or k c 
This is tantamount to changing 



a (l-k)k 2 -(l-k)(k-q) (l.k)q 2 -(l-q)(k-q) 
1 ' ' " 19 fcV - • q) 2 



k 2 q 2 — (k ■ q) 2 
under the integral sign. 

The final result is given by 

= 4e 2 e^ Xa q x k a [A(k, q) + A(k, -q)} 
dH (l-q). V(q) 



(3.5) 



(3.6) 



with 



A(q,k) 
V(q) = 



i(27r) 4 (l-q) 2 (l-k/2) 2 
1 



[B + (k-l)C/2] : 



k 2 q 2 — (k ■ q) : 



[^(k ■ q - 2q l ) + q fl (2k ■ q - k" 



(3.7) 



One can show that for q\ — q\^ q 2 and \q 2 \ — ► oo the term (I — q) 2 can be replaced by q 2 . 
We shall analyze below in detail when it is safe to make this replacement. For the moment 
let us accept it and notice that eq.(|3.6|) then involves 

[V»(q)-V»(-q)}-q» = -2, 

and 



\(q-l)(q-k)-q 2 (k-l)l 



k 2 q 2 — (k ■ q) 2 

and the remaining terms in the integral do not depend on q, so that effectively l x 
k x (k ■ l)/k 2 under the integral sign. 

Taking all this into account we readily get 



Mp V {q, k) = 8e 2 e^ x *q X k a \ 

q z 



dH 



1 



B + \(k-l)C 



1-2 



(k-l) 
k 2 



i(2n) 4 (l-k/2) 2 l~ ' 2 V " " k 2 )' ^ 

We see that up to a numerical factor the last integral is nothing other than in eq. ( |3.2| 
Comparing with eq.(l.l) we see that 

2e 2 f n 



for qf = q 2 = q 2 and \q 2 \ - 
given in eqs.(1.3) and (1.4). 



M(qlq 2 2 )--^-^ 
oo, which coincides with the OPE result 



for q^ 



(3.9) 



-oo 
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4. Asymptotics for q\ ^ q\ 

To analyse the situation at q\ ^ q\ we choose the pion rest frame as a reference frame. 
If we also choose the z axis to be directed along q = (cfi + <?2)/2, one has for the time and 
third components of k^ and q^ 

k» = (k , k 3 ) = (r?v,0), ^ = (go, Qs) = (k • q/m w , ^{k ■ q) 2 - (4-1) 



where 



q 2 = \{ql + ql) - \™l « ^(5? + d). (fc • q) = \(ql - ql)- 



(4.2) 



We are interested in the limit \qf |, \q 2 \ 00 a t fixed u, where 

(k -q) ql- ql 



2 i 2 ' 
91+52 



-1 <u < 1. 



(4.3) 



The condition |w| < 1 follows from our assumption that q\ and g| are of the same sign 
which is the case of physical interest. In this limit one has for k 11 and q^ 



P = m,(l,0,0,0), ^ = i^i (1,0,0,1). 



m 7 



(4.4) 



One sees that, since in our reference frame the components of are of the order of hadronic 
scale (~ 1 GeV), the scalar product (/ • q) is of the order of (k ■ q), or more accurately 



(*-g) = ^(io-i 3 ). 



m 7 



(4.5) 



Using this we can simplify the general expression ( |3.6| ) at (k ■ q) 2 » k 2 q 2 (k 2 = m 2 ) and 
obtain 

M = 4e 2 e^ Xa e^ qi )e"(-q 2 )q x k^x (4.6) 
dH B + (k-l)C/2 



1 1 

(T^F + WW- 

In our reference frame eq. Q4.6| ) can be rewritten in the form 



z(2tt) 4 (/ - k/2f 



1 - 



M 



^e 2 e^ Xa e^q 1 )e"(-q 2 )q x k 17 



2(1 ■ q) 
(k-q) 



(4.7) 



x 



dH B + mJoC/2 
z(2tt) 4 I 2 - mJo 



1 



+ 



1 



1 - 2(Z - l^u/m* 1 + 2(l - Z 3 )u;/m 7r 

2(z -z 3 ; 



x 1 



TOti 
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In the last expression we neglected l 2 /q 2 << 1 in the denominators in the square brackets 
and in the first denominator. 

Consider now the limit m n — ► (the chiral limit). At first glance the amplitude A4 
has a natural dependence on the parameter ui / m n and one might expect a rapid variation 
of the function M. at uj ~ m^/A << 1, where A ~ 1 GeV is a characteristic hadronic 
scale. This is not, however, correct. 

Expanding the denominators in the square brackets in eq. (|4.7| ) we get a series in 
the parameter [u(Iq — l^/m^] 2 . It may seem that the terms of this series are singular at 
— > 0. However the integrals multiplying all negative powers of m n vanish. 

To see this let us expand all the other factors i.e. the denominators (Z 2 — m. K lo)~ 1 ajid 
the functions B and C which depend on (k ■ I) = m n lo, in the parameter (m^lo). The 
integral Q4.7Q becomes 



x 



dH 

z(2tt) 4 



E< 

.fc=0 



M = 8e 2 6^ Xa e»( qi y(-q 2 )q A k°x 

2k oo 



2k ( 2 (^o - k 



-E- 

k=0 



2k ( 2 (^o - h 



2k+l 



X 



p=0 



where we have introduced the expansion 



B + mJoC/2 ^ 2 

L '""7T<"0 n 

p=0 



(4.9) 



We do not specify at the moment the functions A p (l 2 ). In the expression ( f4.8| ) one can 
rotate Iq U4 since the singularities of A p (l 2 ) in Iq obey the Feynmann rules; their 

positions in the complex Iq plane are determined by Iq = ± + 1 2 — ie where A is a 
mass parameter. This leads to the integration in euclidean space 



dH 



E 



dH 



i(27r) 4 (2tt) 4 



d Ie = dl^dl^d l±, d l± = dl±dl2- 



(4.10) 



We now introduce 



U = h \ cos (f), I3 = l\ \ sin . 



(4.11) 



and see that the nonvanishing contributions in the r.h.s. of eq. ( |4.8|) come only from the 
terms with p > 2k in the first sum over p and with p > 2k + 1 in the second one. Indeed 



(4.12) 
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and the integration over <p yields zero contribution unless the power of Iq is larger than or 
equal to that of (Iq — I3). 

Therefore the amplitude can be represented as a double series in uj 2 and m 2 . We shall 
now try to calculate the leading term in m 2 for which purpose one should leave in eq. Q4.8Q 
only the contributions from p = 2k and p = 2k + 1. 

One has 



x 



dl% 



M 



dlt 



x 



2e 2 e MI/A<T e"(g i )e' / (-g2)g A fc g 
4n 2 q 2 

00 00 

J2 u 2 \iif k A 2k + ^ 2fc (/f|) 2fc+1 A 2fc+1 



(4.13) 



,fc=0 



k=0 



The integrals over If are readily calculated and one obtains 



M = 2e%^( qi )e^-q 2 ) q ^ [°° ^ ^ ( (l 2 E ) 2k+1 A 2k + (l 2 E ) 2k + 2 A 2k+1 ^ 2/ „ 



(4.14) 



4ttV J n f-i V 2A; + 1 ' 2A; + 2 

To determine the coefficients A 2 /c and A 2k+ i we expand the factors entering eq. ( ^.9|) 

B ; t!'^ /2 = t^+lmJ C p )( m J^±^, (4.15) 



<J = 



d 2 Y 



where we have used the fact that B and C are the even functions of (k ■ 1) = tu^Iq. The 
coefficients A 2k and A 2k+ i in ( [4.14] ) are equal to 



A 



2k 



E 



Sp 1 C p 



^ (/ 2 )2fc -2 P+ l 2 ^ o (/2)^-2p' 



(4.16) 



A 



2k+l 



E 



1 ^ 

^E 



(l2\2k-2p+2 1 2 ' (72\2fe-2p+l " 



It is implied that in the expression (|4.16| ) for A 2k the terms proportional to C p are absent 
for k = 0. 

Substituting (|4.16| ) into ( f4. 14| ) and changing the order of summation in p and k one 
can derive the following expression for M. 



KA _ 2e 2 e^ A(J e^(Qi)e"(-Q 2 )Q A ^ 
y w — z — ~ — ~ x 



8n 2 q 2 



(4.17) 



x 



E 

p=0 
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/°° / I l 2p 

di% (i 2 E ) 2p+1 c p (i 2 E ) UjL^-Mu) 



where 



^") = E (t+1)( a +1) . (" 8 > 

<t>o{u) = </>(u) = -J [(1 + w)ln(l + w) + (1 - w)ln(l - w)]. 

Generally the coefficient functions B p {l 2 E ) and C p {l 2 E ) in the expansion of S(/ 2 , (& • Z) 2 ) 
and C(/ 2 , (k ■ I) 2 ) in (/c ■ I) 2 are unknown. However we shall see that if the form factors 
B and C decrease fast enough at infinity as functions of the quark masses (virtualities) 
m\ = (/ + k/2) 2 and m 2 = (Z — /c/2) 2 , then the coefficients B p and C p decrease with p and 
since also the (f> p (w) decrease with p one can obtain a reasonable estimate by keeping only 
the first term in the sum in eq. ( f4.17| ). 

To see this we consider a simple model when the mass dependence of B (or C) is 
determined by 

B = ( jir, „2wa9 =5t 1 ( 4 - 19 ) 



(A 2 - m 2 )(A 2 - m 2 ,) 

where we omit an inessential constant factor; and n > 1 is an integer. 
For this model the coefficients B p (l E ) read 



n(n + 1) • • • (n +p — 1) 1 

"pi (a 2 + /|) 2 -+ 2 p' 



- n 7To. , ,2^^^ ' P ^ ( 42 °) 



a 



(A 2 + i|) 2 "' 
The integrals in eq. ( [4.17] ) are of the form 



f oo 

I p = (A 2 ) 2 - 1 / dl% tiB p {l\) = (4.21) 
Jo 

_ n(n + 1) • • • (n + p- 1) (2p)! 

~ pi (2n - l)2n(2n + 1) ■ • ■ (2n + 2p - 1) ' 

Therefore 

11 3 

7 ° ~ 2n - 1' 7l ~ (2n- l)(2n + 1) ' /2 ~ (2n - l)(2n + l)(2n + 3) ' < ^' 22 ' > 

15 



(2ra- l)(2n+l)(2n + 3)(2n + 5)' 

We see that the ratio ii//o = l/(2n + 1) is only 30% for n = 1 and 20% for n = 2, while 
the other coefficients are even smaller. 
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Keeping only the p = terms in eq. Q4.17 ) one obtains 



M = — 2 [f n (f)(u)/3 + g^l - <p[u))\, (4.23) 

where 

1 r°° 

9* = W J ^eC^e) (424) 
and we have used the fact that the expression ( |3.2| ) for simplifies in the chiral limit to 

u= ^r + i i z c °^) ■ (4 - 25) 

Eq. (|4.23|) provides an expression for the amplitude M. in terms of one unknown parameter 
g n and, as such, is of some interest. We shall argue however that g n « and thus the 
term proportional to gf w may be neglected. In the next section we show that -B(O) = 
1/ f n . Thus -B(O) describes the low energy interaction of pions with quarks (see eq. (1.9)) 
reflecting the Goldstone character of pions. The form factor C(0) does not show up at all 
in this interaction at low energies if the PCAC ideology is used. Since = 93 MeV is 
"abnormally" small in the units of usual hadronic scale A ~ 1 GeV, the value of B(0) is 
"abnormally" large, but we would expect that C(0) has a "natural" value C(0) ~ 1/A 3 . If 
this is correct then the term involving C in eq. Q4.25|) should be negligible and we would 
expect roughly 

A 2 i?(0) _ A 2 



(4.27) 



Then for g n in (|4.24|) we expect 

A 4 C(0) J^^U . 
9lT ~ 8tt 2 ~ 8tt 2 A ' U 

so that f n ~ lOg-v- 

Neglecting the g^-term in eq. Q4.23| ) we arrive at the final result 

1 00 2k 

m = - 2 \a + W )in ( i + u) + (i - »w - «)] = E (t+1 ) (afc+1) 

which was quoted in Introduction. 

In getting eq. ( |4.28|) we neglected the second term ~ g n and all terms with p ^ 
of eq. ( f!.17| ). Since for n = 2, which seems a reasonable value in the model (4.19), the 



accuracy of the second of these approximations is about 20% and there is an additional 
smallness related to the fact that <p p+ i(l) < O.30 p (l) we expect the accuracy of the final 
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result ( fl28| ) might be of order < 10%. The value </>(0) = 1 corresponds to the OPE result 
for the asymptotics at q\ = q\. On the other hand, 0(1) = 2 In 2 « 1.4. This has 
to be compared with the result 0(1) = 3/2 given in ref. (see Introduction). Since we 
used certain numerical approximations and our whole approach is different from that of 
ref. the 7% coincidence of these two figures seems to be quite satisfactory. We plot 
this function 4>(uj) in Fig. 10. It is interesting to compare our result for the function 0(w) 
with that given in ref. ||. In notation of ref. the function (denoted there as I{oj)) is 
represented by an integral 



In the limit of infinitely large renormalization scale /i^oowe have | || | [£L4 



^(^ 2 )-0T(O = ~(i-£ 2 )- 

As it is pointed out in ref. || this expression for </u(£, fj 2 ) is reliable only for \u>\ < 1/2. To 
illustrate the difference between the functions and I we plotted the ratio 2(0 — I)/(<p + 1) 
in Fig. 11. It is easy to see that for \u\ < 1/2 the difference is 1% . At \u\ = 1 the 
difference is 8% which is not surprising since we use quite different physical assumptions. 
Thus the above estimate of our accuracy is compatible with the difference of our result 
with that of ref. 131 . 



5. 7T — > 27 and the anomaly 

It is interesting to understand how the formalism developed in this paper corresponds 
to the description of n° decay into real photons by PCAC with the Adler-Bell-Jackiw 
anomalous term. We shall see below that the comparison of these two approaches leads to 
the relation 

B(0) = (5.1) 

where B(0) is the value of the form factor B with all momenta zero. 

First of all it should be stressed that since q 2 = (qi + g 2 ) 2 /4 = —m 2 /4 is now small 
there is no reason to consider only the triangle diagram which has been used above. The 
most general amplitude of n° — > 27 decay can be represented as shown in Fig. 8. The 
A^2 amplitude of Fig. 86 is of the most general form provided that the four-point Green 
function v v is of the most general character. In the M.2 amplitude we have arbitrary 
chosen to show explicitly the quark which emits the (q±, fj) photon while the vertex for the 
emission of the second photon, (52, v), is not shown and is absorbed into v v . Separating 
the Ai\ amplitude (Fig. 8a) from A^2 means that we exclude the pole contribution to v v 
depicted in Fig. 9. The reason for doing so is that in this section we shall intensively use the 
chiral limit, together with Adler's condition for the emission of soft pions. This condition 
requires that any amplitude should vanish in the limit when the pion four-momentum goes 
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to zero except the pole terms which can survive. Because of Adler's condition v v being 
defined without a pole contribution vanishes at /c M — > 0, i.e. we have 

v v -»• at fc M -> 0. (5.2) 

The idea of the argument below is to prove that in computing the amplitude Ai 2 we can 
use the Green function v v at q 2 = 0. Given this we can use the Ward identity for the 
emission of a photon with zero momentum 

^(M,0) = -^. (5.3) 

Thus we shall manage to compute A4 2 in terms of the T(k,l) vertex used throughout this 
paper. 

We now give to the proof that v u (k, I, q 2 ) can indeed be taken at q 2 = 0. 
The most general structure of v v which survives in the calculation of the amplitude 
M 2 is 

v v (k,l,q 2 ) = (F lu k , + F 2 J + F 3u fi 2 +F 4 <y u )<y 5 , (5.4) 

where all functions Fi„, F 2v , F% v and F4 depend on the vectors k, Z, and q 2 . 

Calculating the trace appearing in computation of the amplitude A4 2 (Fig. 86) one 
obtains the following expression 

e xp ^[F lv l p q 2(T k x + F 2u k p q 2 J x /2 + F 3 J p k a q 2X + F 4 (l p k a - l p q 2a + k p q 2a /2)]. (5.5) 

Consider each of the terms in ( |5 . 5| ) . 

For the F\ v term we already have the necessary structure to provide the e PiV \ a e p, {qi) 
e u (—q 2 )q x k (T factor (since e\ PPiC7 q 2cr k\ = e\ p ^ a q 2cr qi\). Therefore the contribution propor- 
tional to F\ v is already of the second order in the small momenta (photon momenta) and 
thus all the other factors can be taken at zero momenta. Particularly we can take F\ v at 

g 2 = o. 

For the F 2v term the situation is basically the same and again F 2v can be taken at 
q 2 = 0. In this case, however, F 2v does not actually contribute since in Adler's soft pion 
limit F 2v — > as k^ — > 0. This follows from an application of Adler's soft pion theorem to 
F 2u {k,l,q 2 ). 

The F^v form factor can also be taken at q 2 = but now this is not the same as 
putting q 2 = in v v (unlike the case of the first two terms in ( |5.4| )) since F 3t/ (/c, /, 0)^2 
in principle could survive. Luckily the whole contribution to M. 2 , proportional to F 3i ,, is 
again of the third order in the small momenta. Indeed k a q 2 \ — > q\ a q 2 \ and also F 3l/ is 
itself proportional to the pion 4-momentum due to the soft pion emission theorem. 

Lastly the F 4 form factor can also be taken at q 2 = 0. Indeed part of it, proportional 
to (k-q 2 ) is immediately negligible, while the part oc k-l can contribute through F 4 l p q 2(J — > 
F4qi p q 2 a, but then F 4 can obviously be taken at q 2 = 0. 
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Thus we see that for different reasons all the contributions to v v in eq. (|5.4| ) can be 
taken at q 2 = 0. Therefore v v (k,l,q 2 ) — > v u (k,l,0) and the fundamental Ward identity 
( p73|) is applicable to our calculation. 

Using for T(k, I) the expression 



F(k,l) = [iBJk+i(k-l)C}] l5 

and differentiating it w.r.t. l v we need to keep only the terms 

dB{l 2 ^ 



dl< 



+ (k ■ /)C 7 ,]75 



(5.6) 



(5.7) 



as follows from the consideration given above: only the F\ v and F± form factors can actually 
contribute to the amplitude in question. 

We are now ready to write down the final result for A4i and Ai 2 of Fig. 8 which can 
be obtained by a straightforward calculation. For A^2 we obtain 



A 1,0 



2e 2 e^ Xa e^(q 1 )e-(-q 2 )q A k 
8tt 2 



m - 1 / Micvi) 



(5.8) 



For the amplitude A^i the result is 



Mi 



2e 2 e^x^( qi )e^-q2)q A k 
16tt 2 



dl%C(l%). 



(5.9) 



Thus the sum of M.\ and A4 2 does not contain C and there results 

2e 2 e^ Xa e^ qi y(-q 2 )q x k° 



M= Mi+ M 2 



8tt 2 



B(0). 



(5.10) 



We can compare this result with the usual PC AC amplitude 

2e 2 e^ Xa e»{q x y{-q 2 )q x k° 



M 



to get the relation 



8tt 2 A 



1 



B(0) = T . 

Jit 



(5.11) 



(5.12) 



This equation can be understood in the following way. In the chiral limit of vanishing 4- 
momentum of the pion the pion-quark-antiquark interaction has a pseudovector character 
(corresponding to the 5-form factor) with a fixed coupling constant 



L e ff = -r{d^^)[u^^u - d-y^d}. 

J7T 



(5.13) 



Equation (|5.13| ) reveals the Goldstone character of the pion. 
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6. Conclusions 



We have studied the 7r°77 vertex for a range of virtual photon squared masses qj,q%, 
in a rather general picture where the pion is treated as a composite, non-pointlike particle. 
We show how our approach reproduces the operator product expansion results valid for 
Qi = <?2 ~ * 00 an d compare with the infinite momentum frame result for q\ = 0, q\ — > 00. 

We have been able to find a closed analytic formula valid for arbitrary fixed values 
of the ratio q\jq\ in the limit that at least one of \q\\ and \q\\ is large. If the Feynman 
amplitude is written 

M = M(qlqi)e^ Xa e^q 1 )e"(-q 2 )q x ^ 



where q = (qi + <?2)/2, k = q\ — qi and if we define 

2 2 

q{ - qj 
ql + ql 

then our result is 



2 2 

UJ = s-, — 1 < < 1 



with 



M(^ 2 2 )=M(O,O)-(-^0^) 0(u/) 

= -^[(1 + w) log(l + log(l - a;)]. 

The amplitude is of great interest theoretically because of its link to the axial anomaly. 
Indeed this connection has led to dramatically incorrect estimates of the amplitude at large 
photon virtualities. 

On the experimental side the amplitude M plays a crucial role in several physical 
processes. For q\ < 0, q\ < it is relevant to the two-photon physics reaction e + e~ — > 
7r°e + e~, and for q\ > 0, q\ > it controls the annihilation process e + e~ — > 7r°/ + /~. It 
would be of great interest to test our prediction experimantally. 
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Figure captions 

Fig. 1. Feynman diagram for e~ + e~(e + ) — > e~ + e~(e + ) + n° reaction. 

Fig. 2. Feynman diagram for e~ + e + — > e~ + e + + n° and e~ + e + — > \i~ + + n° 
reactions. 

Fig. 3. Proper vertex for 7r° — > qq decay. 
Fig. 4. Diagram for calculation of f^. 

Fig. 5. Diagram for calculation of the n° — > 7* + 7* amplitude. 

Fig. 6. Possible gluon corrections to the diagram of Fig. 5. 

Fig. 7. Possible radiative corrections to the diagram of Fig. 5. 

Fig. 8. a), b) Diagrams for contributions to the tt° — > 27 amplitude. 

Fig. 9. Diagram for the pole contribution to the amplitude. 

Fig. 10. Plot of the function 4>(u). 

Fig. 11. Plot of the ratio A(u) = 2(<f>(u) - I(u))/(<j>(u) + where I(u) is the 

result of ref . . 
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Fig.8 

27 




28 



0.2 0.4 0.6 0.8 1 

Fig.10 



29 




30 



